Effects of the anomaly on the two-flavor QCD chiral phase transition 
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We use strongly coupled lattice QED with two flavors of massless staggered fermions to model 
the physics of pions in two-flavor massless QCD. Our model has the right chiral symmetries and can 
be studied efficiently with cluster algorithms. In particular we can tune the strength of the QCD 
anomaly and thus study its effects on the chiral phase transition. Our study confirms the widely 
accepted view point that the chiral phase transition is first order in the absence of the anomaly. 
Turning on the anomaly weakens the transition and turns it second order at a critical anomaly 
strength. The anomaly strength at the tricritical point is characterized using r = {M^i — M-jr)/ Pn' 
where M^/ , M,r are the screening masses of the anomalous and regular pions and p^' is the mass-scale 
that governs the low energy fluctuations of the anomalous symmetry. We estimate that r --^ 7 in 
our model. This suggests that a strong anomaly at the two-flavor QCD chiral phase transition is 
necessary to wash out the first order transition. 
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Understanding the nature of the finite temperature chi- 
ral phase transition in QCD has been an area of active 
research for almost two decades. The subject is impor- 
tant for the field of relativistic heavy ion experiments [ij . 
It also plays a role in cosmology where a first order transi- 
tion can lead to rich cosmological scenarios Q ■ The con- 
ventional wisdom regarding the nature of the chiral phase 
transition in A^^-flavor QCD comes from the renormaliza- 
tion group study of linear sigma models which describe 
the fluctuations of the order parameter in QCD. These 
models are invariant under SU{Nf) x SU{Nf) x U{1) chi- 
ral symmetry. The anomaly is introduced through a term 
that breaks the C/(l) symmetry explicitly. These studies, 
usually based on the e-expansion reveal that the number 
of light quark flavors and the strength of the anomaly 
can play a significant role in determining the order of the 
transition [3|. While a second order phase transition is 
possible for Nf = 2, the presence of more light quarks 
can introduce fluctuations that will force the transition 
to become first order. Similarly, if the anomaly is suffi- 
ciently weak at the transition, the two-flavor transition 
could itself turn into a first order one. It is known that 
the anomaly weakens with increase in temperature and 
number of colors. What really occurs in nature remains 
a topic of active research even today. [1, d, Recent 
reviews can be found in 0, Q • 

Although the e-expansion is known to be unreliable in 
three dimensions, Monte Carlo methods have confirmed 
some of its predictions. Both non-perturbative studies 
in the sigma models and lattice QCD calculations have 
found that for three or more flavors of light quarks the 
chiral transition is flrst order as predicted. Further, there 
is little doubt that in the presence of a large anomaly, the 
two flavor transition can be second order in the three di- 
mensional SU{2) X SU{2) ~ 0(4) universality class 0. 
On the other hand, the effects of a weak anomaly re- 
mains relatively unexplored. Studies within mean field 
theory, not surprisingly, agree with the predictions of the 



e-expansion, which forbids the two-flavor transition with- 
in, Hi- However, 



out the anomaly to be second order 
recent renormalization group studies, based on three di- 
mensional perturbation theory and resummation tech- 
niques, suggest that a second order phase transition is 
indeed possible in the absence of the anomaly [13] ■ It 
would be interesting to identify and study this new, 
SU{2) X SUi2) X f/(l) ~ 0(4) X 0(2) second order critical 
behavior using Monte-Carlo methods if it exists. To do 
this, one will have to study 0(4) x 0{2) symmetric sigma 
models with the relevant symmetry breaking pattern and 
look for the predicted continuous transition. As far as we 
know this is difficult in the conventional formulation of 
the sigma models due to the lack of efficient algorithms 
and hence remains unexplored. 

Here we offer a new approach to the subject. Using 
strongly coupled lattice QED with two ffavors of stag- 
gered fermions, which can be studied ver y ef ficiently with 
recently discovered cluster algorithms [l3|, ll5|, we 
model two-flavor massless QCD in the absence of the 
anomaly. Our model has the right symmetries and sym- 
metry breaking pattern and so should be viewed as an 
alternative formulation of the relevant 0(4) x 0(2) sigma 
model. An additional four-fermion term can be used to 
introduce the effects of the anomaly and this allows us 
to efficiently study the phase diagram as a function of 
temperature and anomaly strength using Monte Carlo 
methods for the first time. Although our study is not an 
exhaustive search for the second order critical point pre- 
dicted in 1^, it can be viewed as a first attempt within 
the context of a specific model. 

The action of our model is given by 



S 



X l_L—l 



(1) 
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where x = {x, t) denotes a lattice site on a L'^ x 4 hyper- 
cubic lattice. Va; are two component Grassmann 



nature of the fermion fields the following constraints must 
also be satisfied at each site x: J^^'^lii^) — 0' 2/(x) -|- 
fields that represent {u,d) quarks, and 0^,^ is the com- T,JK(^) + '^i(^) + \K(^)\^ = 2, T, JKi^) - '^ti^)] = 0> 



pact U{1) gauge field through which the quarks inter- 
act. Here /i = 4 denotes the temperature direction, while 
fi — 1,2,3 are the spatial directions. The staggered 
fermion phase factors 77^, a; obey the relations: rj\ ^ — T 
and 77^^ = 1 for i = 1, 2, 3. The parameter T controls the 
fictitious temperature. The coupling C sets the strength 
of the anomaly. 

The above action exhibits a global SUl{2) x SUr{2) 
symmetry. Indeed it is easy to check that the action 
is invariant under ijjx^ —>■ Lipx^^i^xo ^ R'^Xo,^x ~* 
^x„L^T^x^ ~^ i'x^R^ 7 where L,R G SU{2) and Xe and 
Xo refer to even and odd sites. When C — the action 
is also invariant under a C/(l) transformation given by 

i^x^ e'^i^x,:i^x^ &~'^-<Px^,'^xa '>Px^e-''^,'^x, 
■ip,j.^e^^ . This will be interpreted as the anomalous /7a(1) 
transformation of QCD. When C ^ 0, the J7a(1) is ex- 
plicitly broken to Z2 and the action is invariant only un- 
der SUl[2) X SUr{2) X Z2. Thus, our model has the 
same symmetries as — 2 QCD with a parameter C 
which helps change the anomaly strength. 




FIG. 1: An example of a DPI configuration in two dimensions. 

The partition function of our model can be expressed 
as a classical statistical mechanics problem involving 
gauge invariant objects: dimers (D), pion loops (P) 



Djects: 

and instantons (1)0, [H, ll7|. We refer to the set of 
these configurations as DPI configurations and denote 
them {/C}. Each DPI configuration is characterized by a 
site variable I{x) S {0,1} representing instantons and 
three bond variables, 7r^(a;) G {0,1}, 7r^(a;) G {0,1}, 
{^170,1}, which give rise to dimers and pion 
= 7r^(a; — ji) and simi- 
Due to the Grassmann 



loops. In our notation 7r"^(a:) 
larly for other bond variables. 



±4. 



where the sum over /i goes over ±1 
an illustration of a DPI configuration in 1 



Figure [T] gives 
- 1 dimensions. 



The partition function is given by Z = ^j^j e ^ where 
S = + <{x) + 7ri(x)]log(r) + /(x)log(C). 

Although the SU{2) x SU{2) x Ua{1) symmetry of the 
model is not apparent in the DPI formulation, it is easy to 
check that Jlix) = ^i(x), J^(x) = e(x)[Ti^^{x) - 7r^(a;)] 
and J^{x) = £{x)[K'^{x) + n{{x) + \nl{x)\-ll2\ form the 
vector, chiral and axial conserved currents when C — 0. 
Note £{x) is defined to be -fl on even sites and —1 on 
odd sites. 

We will focus on four observables: 



(1) Vector current susceptibility 
3 / 



c 



3Z^' 



1=1 




(2) 



(2) Axial current susceptibility 



(3) Chiral condensate susceptibility 



(3) 



(4) 



(4) Mass of the anomalous pion M^/ obtained from the 
zero momentum correlation function 

^ ^ ^ {e^-ip^x^XA^) eip'ilj{Q,y\_)), (5) 

by fitting it to the form Acosh(Af^x). Here x± 
and 2/_L represent coordinates perpendicular to the 
direction x. 

At C = and small T we expect that the SUl{2) x 
SUl{2) X C/a(1) symmetry of the model is broken to 
a diagonal SU{2) flavor symmetry. As T is increased 
one expects a phase transition at T = Tc where sym- 
metries are restored. This phase transition can be stud- 
ied conveniently using the two spin-stiffnesses defined as 
= hmL^oo Yc and p,,' = limi^oo Ya- We expect 
these two quantities to be non-zero in the broken phase 
and to vanish in the symmetric phase. With knowledge 
of these, one can write down the leading term of the low 
energy effective action that governs the Goldstone-boson 
fluctuations in the broken phase. This action turns out 
to be 
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(fa 



P-n 



Tr 



(6) 
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FIG. 3: Plot of LYc versus T for different values of L for C = 
FIG. 2: Plot of LYa and LYc versus T for different values 0.3. The presence of a point where aU the curves cross shows 
of L for C = 0. The lack of a point where all the curves the presence of a second order transition. We can estimate 
cross shows the absence of a second order transition. We can = 2.83555(10). The inset shows that all our data fits 0(4) 

estimate Tc ~ 2.466(1) scahng extremely well. 



with U{x) e SU{2) and u{x) — exp(iri'). Clearly and 
Pri' define mass-scales for the problem. If the transition 
to the symmetric phase is second order, we expect p ~ 
A{Tc—T)'^. We also expect the finite size scaling relation: 



LY{L) ^ ^ /, [(T, - T)ii/-' 



k=0 



(7) 



in the region (Tc — T)V-/'^ <^ 1. In other words a second 
order transition will be characterized by the fact that 
LY{L) is independent of the lattice size at T = Tc- We 
will use this characteristic to decide if a transition is 
second order or not. In Figure [2] we plot LYa{L) and 
LYc{L) as a function of T at C = 0. Although there 
is clear evidence for a transition at Tc ~ 2.466(1), it is 
not second order. This result confirms the e-expansion 
scenario. 

As discussed in the introduction, when the anomaly 
is large one expects 0(4) critical behavior. In order to 
test this we set C — 0.3. In this case we expect crit- 
ical behavior only in Yc{L) since the anomalous pion 
irj\ will remain massive. In Figure [3] we plot LYc{L) 
as a function of T. Unlike the C — case, now all 
the curves cross at a point. We can fit all our data 
to the 3d 0(4) second order scaling form {v = 0.745 
0) given in eq.d?]) with a x^/DOF = 1.3. The fit 
yields /o = 0.494(2), h = 0.388(4), ^ = 0.022(4), 
/3 = -0.008(2), Tc = 2.83555(10). The scaled data is 
shown in the inset of Figure [H 

Given that the transition is first order at C = and 
second order at C = 0.3 it is interesting to locate the 
tricritical point. After some trial and error we have found 
that the tricritical point is located at C « 0.030(5) and 
Tc ~ 2.52652(3). At a tricritical point one expects mean 
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FIG. 4: Plot of LYc (top) and x/L^ (bottom) versus T for 
different values of L for C — 0.03. The scaling of x assumes 
77 = 0. The presence of a point where all the curves cross 
shows the presence of a second order transition. We can esti- 
mate Tc ~ 2.52652(3). 



field scaling in three dimensions [l8|, i,e.. v — 0.5 and 
77 = 0. In addition to eq. here we have also looked 
at the critical scaling relation for the chiral condensate 
susceptibility 



L'-^Y.9k[iTc-T)L 



(8) 



fc=0 



In Figure [4] we show the behavior of LYc and x/i^ as a 
function of T for C = 0.03. Clearly, again the curves meet 
at a point. A combined fit of all our data to eqs. (|7l8p . 
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FIG. 5: The plots shown are close to the tricritical point 
C = 0.03 and T = 2.5265. Assuming u = 0.5 and r? = 0, 
in the top panel we show LYc and x/L^~^ ^ function of 
(T - Tc)V-^" for all of our data. The fact that all the data fall 
on a single curve demonstrates mean field universality. The 
bottom left panel is the plot of Ya versus L and the right 
panel shows the correlation function G{x). Solid lines are fits 
discussed in the text. 



assuming mean field scaling, gives a /DOF = 1.1. The 
fit yields /o = 0.853(3), /i = 0.091(2), go = 0.384(1), 
gi = 0.078(2), Tc = 2.52652(3). In the fit, fk^gt for 
k > 2 were set to zero. In the top plot of Figure [5l we 
show the scaled observables. 

It is useful to characterize the strength of the anomaly 
along the critical line from C = to C = 0.03 using a 
dimensionless parameter. In order to do this we consider 
two mass scales naturally defined along this line. The 
first mass scale is pr,' = [(PjjOt^ + (/'»;')t+]/'^ which is 
the average value of the anomalous spin stiffness across 
the transition. The second mass scale is the difference 
between the screening masses of the anomalous pion 
Mri' and the regular pion M^. The dimensionless ra- 
tio r = {Mfji — Mt^)/ Pj^i clearly measures the strength of 
the anomaly: It is zero at C=0, i.e., in the absence of the 
anomaly, and increases with C. In Figure[5] (bottom- left) 
we show the finite size scaling of Ya{L) close to the tri- 
critical point. By fitting the data to the form p'^ + a/ L we 
find that p'^ w 0.05. Note {pri')T- = iPv')T^ ^^"^ 
critical point. We have noticed that p^' does not change 
much from C = to C = 0.03 along the critical line. 
This suggests that r is indeed a good parameter to mea- 
sure the strength of the anomaly. In Figure [5] (bottom- 
right) we show the correlation function G{x) again close 
to the tricritical point, for L = 20. By fitting the data 
to Acosh(M^a;) we estimate « 0.35. Since M^r = 
at the tricitical point, we conclude that r w 7. Thus, 
r changes from zero in the absence of the anomaly to 
approximately 7 at the tricritical point. 



Although r is an interesting ratio, it is clearly non- 
universal since it involves non-critical scales. So we can- 
not predict the value of r in QCD from our results. How- 
ever, the large value found here at the tricritical point 
may be a generic feature of many theories including QCD. 
Thus, a strong anomaly may be necessary before the sec- 
ond order 0(4) critical behavior sets in. On the other 
hand lattice QCD calculations suggest that the anomaly 
is weak at the chiral phase transition [l3| • If this is true 
a second order chiral transition in the 0(4) universality 
class for two-flavors of massless quarks in QCD seems 
unlikely. 
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